Mott-Hubbard exciton in the optical conductivity of YTiOa and SmTiOs 
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In the Mott-Hubbard insulators YTiC>3 and SmTiOs we study optical excitations from the lower 
to the upper Hubbard band, l^d 1 ) — > \d°d 2 ). The multi-peak structure observed in the optical 
conductivity reflects the multiplet structure of the upper Hubbard band in a multi-orbital system. 
Absorption bands at 2.55 and 4.15 eV in the ferromagnet YTiOs correspond to final states with a 
triplet d 2 configuration, whereas a peak at 3.7 eV in the antiferromagnet SmTiC>3 is attributed to 
a singlet d 2 final state. A strongly temperature-dependent peak at 1.95 eV in YTiC>3 and 1.8 eV 
in SmTiOs is interpreted in terms of a Hubbard exciton, i.e., a charge-neutral (quasi-)bound state 
of a hole in the lower Hubbard band and a double occupancy in the upper one. The binding to 
such a Hubbard exciton may arise both due to Coulomb attraction between nearest-neighbor sites 
and due to a lowering of the kinetic energy in a system with magnetic and/or orbital correlations. 
Furthermore, we observe anomalies of the spectral weight in the vicinity of the magnetic ordering 
transitions, both in YTi03 and SmTi03. In the G-type antiferromagnet SmTiOs, the sign of the 
change of the spectral weight at Tjv depends on the polarization. This demonstrates that the 
temperature dependence of the spectral weight is not dominated by the spin-spin correlations, but 
rather reflects small changes of the orbital occupation. 

PACS numbers: 71.20.Be, 71.27. +a, 77.84.Dy, 78.20.-e, 78.20.Ci 



I. INTRODUCTION 

In strongly correlated electron systems, the competi- 
tion between kinetic energy and Coulomb repulsion gives 
rise to a variety of intriguing phenomena^ 2 - The most 
simple approach is the single-band Hubbard model, with 
on-site Coulomb repulsion U and a kinetic part propor- 
tional to the inter-site hopping amplitude t. For U larger 
than the band width, the band splits into a lower and an 
upper Hubbard band (LHB and UHB, sec inset of Fig. 
[T]). At half filling one finds a Mott-Hubbard insulator 
with one localized electron per site, i.e., the Coulomb 
energy dominates. However, the low-energy physics is 
determined by the kinetic energy: virtual hopping of the 
electrons to neighboring sites is effectively described by 
exchange interactions with J oc t 2 /U. These govern the 
spin degrees of freedom and, in a multi-orbital model, are 
also relevant for the orbital degrees of freedom4>« 

The competition between Coulomb energy and kinetic 
energy also governs the formation of bound states, e.g., 
excitons. In simple band insulators, binding of an elec- 
tron in the conduction band and a hole in the valence 
band reduces the Coulomb energy, while the kinetic en- 
ergy increases. In Mott-Hubbard insulators, the lowest 
optical "interband" excitation creates an empty site and 
a doubly occupied site, i.e., a hole in the LHB and a par- 
ticle in the UHB. A Hubbard exciton can be regarded 
as a bound state of an empty site and a doubly occu- 
pied site, moving in a background of singly occupied 
sites. Studies of excitons in correlated electron systems 
thus far have focused on one- (ID) or two-dimensional 
(2D) systems. Remarkably, it has been found that ex- 



citon binding can be driven by either the Coulomb en- 
ergy or the kinetic energy. The former is found in the 
ID extended Hubbard model, which takes into account 
the Coulomb interaction V between nearest or next- 
nearest neighbor s it e s£L2uLi£iiLi2i£2ii4 Since both the 
Mott-Hubbard gap and the attractive interaction for ex- 
citon binding result from Coulomb interactions, one ex- 
pects different physics compared to band insulators. In 
fact, excitons are only formed below the gap if V ex- 
ceeds a critical value.— iLSiii For smaller values of V, an 
excitonic resonance is found in the continuum above the 
gap, strongly affecting the line shape of the optical con- 
ductivity o^u;) ^ 11 ' 14 An exciton below the gap has been 
observed in ID Ni-halogen chainS ) 15 i 16 and this exciton 
contributes to the gigantic non-linear optical response 
observed in these compounds. 14 ^ 5 ' 16 ] 17 

The kinetic energy is of prime importance for excitons 
in the 2D cuprates^ i 2 0^., 2 2 1 23,24 1 25 1 26 1 27 1 28 1 29 which arc 

of charge-transfer type. The dispersion of a spinless 
charge-transfer exciton is of order t, larger than the 
single-particle dispersion, which is suppressed to ~ J 
by antiferromagnetic (AF) correlations. Thus exciton 
formation reduces the kinetic energy ) 18 i 19 i 20 i 21 i 22 which 
bears resemblance to a possible mechanism for Cooper 
pair formation in high-T c superconductorsi 21 i 30 i 31 Exper- 
imentally, the exciton dispersion has been studied by 
electron-energy-loss spectroscopy^ and by resonant in- 
elastic x-ray scattering (RIXS) j 28 ' 29 It has been claimed 
that the dispersion is indeed large ) 19 i 28 but recent high- 
resolution RIXS datajSS indicate that the exciton disper- 
sion is suppressed by the coupling to phonons. 

Here, we report on the observation of an excitonic res- 
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FIG. 1: (Color online) Optical conductivity of YTi0 3 at 15 K. 
Inset: sketch of the optical excitations from the lower Hub- 
bard band (LHB) and the oxygen 2p band into the upper 
Hubbard band (UHB) in case of a single, half-filled orbital at 
the transition-metal site. 



onance in the optical conductivity cr(ui) of the 3D Mott- 
Hubbard insulators YTi03 and SmTi03. The former is 
ferromagnetic below T c = 27 K, the latter antifcrromag- 
netic below TV = 53 K, both exhibit orbital order 
Due to the orbital multiplicity in these d 1 spin 5=1/2 
compounds, the upper Hubbard band consists of a series 
of different d 2 multiplets. In YTiC>3, the lowest multiplet 
is identified with a peak at 2.55 eV, whereas a strongly 
temperature-dependent peak at 1.95 eV is attributed to 
an excitonic resonance. For a proper determination of U 
it is essential to take excitonic effects into account. Wc 
discuss the possible relevance of the kinetic energy for ex- 
citon formation in orbitally ordered compounds, similar 
to the case of a 2D antiferromagnet. Our results provide 
the experimental basis to disentangle the role of Coulomb 
and kinetic energy in 3D Mott-Hubbard insulators. 

The spectral weight of the LHB-UHB excitation is 
expected to depend on the nearest-neighbor spin-spin 
correlations i 35 ' 36 i 37 ' 38 In a single-band model, the spec- 
tral weight vanishes in the case of ferromagnetic or- 
der due to the Pauli principle, and one expects a 
strong change of the spectral weight as a function of 
the temperature T at the magnetic ordering transi- 
tion. In a multi-orbital system, the spectral weight 
also depends on the orbital occupation. This kind 
of analysis has been applied to a number of com- 
pounds with different transition-metal ions (Mn, V, 
Ru, Mo)^^SdMl^^i^^ For instance in 
LaMnC>3 and LaSrMn04, a quantitative description of 
the experimentally observed T dependence of the spec- 
tral weight has been obtained i 35 i 38 In the manganites, 
the T dependence is entirely ascribed to the spin-spin 
correlations, whereas the orbital occupation is assumed 
to be independent of T. This reflects the large ligand- 
field splitting A eg of roughly 1 eV of the e g orbitals in 
these d 4 compounds . 35 ! 38 Here, we show that the T de- 
pendence of the spectral weight of YTiC>3 and SmTiO.3 
is not dominated by the spin-spin correlations. This is 
particularly evident for SmTi03 , where the sign of the T 



dependence of the spectral weight depends on the polar- 
ization. This behavior can be attributed to small changes 
of the orbital occupation in these ti g compounds. 

The paper is organized as follows. Section II addresses 
the experimental details. The optical conductivity of 
YTi03 and SmTiOs is reported in section III. In sec- 
tion III. A we first discuss the multiplet assignment and 
argue that the lowest peak has to be interpreted as an 
excitonic resonance in both compounds. A possible con- 
tribution of the kinetic energy to exciton binding in the 
case of antifcrro- orbital order is proposed in section III.B. 
In section III.C wc discuss the temperature dependence 
of the spectral weight and the relevance of spin-spin cor- 
relations and orbital occupation. The anisotropy of the 
spectral weight of the lowest multiplet in YTiC>3 is ad- 
dressed in section III.D. A summary and conclusions are 
given in section IV. The role of oxygen defects for the 
analysis of ellipsomctric data of YTi03 is discussed in 
the appendix. 



II. EXPERIMENTAL 

Single crystals of YTi03 and SmTiOs were grown us- 
ing the floating-zone technique. The crystal quality and 
stoichiomctry were checked by x-ray diffraction, EDX, 
and polarization microscopy. The crystals are single 
phase and single domain. From magnetization measure- 
ments (SQUID, PPMS) we find that YTi0 3 becomes fer- 
romagnetic below T C =27K and SmTiOs antiferromag- 
netic below T/v=53K. Further details on crystal prepa- 
ration and characterization can be found in Ref. 1471 . In 
YTiC>3, four-sublattice orbital order has been reported 
up to room tcmpcraturei 32 i 33 In both compounds an 
orbital-ordering transition has not been observed, i.e., 
they are considered to be orbitally ordered up to the 
melting temperature, or, in other words, the distortions 
arising from the orbital occupation do not break the crys- 
tal symmetry. 



Generalized ellipsometric data— was obtained us- 
ing a rotating-analyzcr cllipsometer (Woollam VASE) 
equipped with a rctarder between polarizer and sample. 
The angle of incidence was 70°. Immediately after pol- 
ishing, the sample was kept in an UHV cryostat. The 
measurement background pressure of p < 1CP 9 mbar has 
been achieved by a bakeout at 400 K for 24 h. Window 
effects have been corrected using a standard Si wafer. 
In orthorhombic i?Ti03, only the diagonal elements a a , 
a b and er c of the complex optical conductivity tensor 
o~(lo) = o"i + i<T2 are finite. In YTi03, we have deter- 
mined er(w) from the normalized Muller matrix elements 
m\ 2 , rn l 21 , m 33 , and TO34, where i = 1 - 4 denotes differ- 
ent orientations of the sample, namely with s-polarized 
light parallel to the crystallographic a and b (a* and c) 
axes on the ab (a*c) surface, where a* = [110] within the 
Pbnm space group. In SmTiOs, (t(oj) has been deter- 
mined from measurements on be and ab surfaces. 

Ellipsomctry is a surface sensitive technique, thus one 
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FIG. 2: (Color online) Optical conductivity of YTiOs in the 
vicinity of the onset of excitations across the gap. Good agree- 
ment is observed between our ellipsometry data and results 
determined from the combination of transmittance and re- 
flectance measurements.— 

has to consider the possible contribution of surface con- 
taminations or adsorbate layers. To this end we have 
polished and measured a sample of YTi03 several times, 
both in UHV and under ambient conditions, and for dif- 
ferent angles of incidence. The raw data show small vari- 
ations which are attributed to the surface. A consistent 
description of all data sets for the two distinct surface ori- 
entations has been achieved by assuming a non-absorbing 
cover layer, where only the thickness d < 2nm of this 
layer has been allowed to vary for different data sets. 
For an extensive discussion of the data analysis, we re- 
fer to Ref. H^. The particular choice of the cover layer 
has a certain influence on the absolute value of o{ui), 
but we emphasize that the temperature dependence is 
hardly affected. We also have checked carefully that the 
observed temperature dependence reflects the properties 
of YTi0 3 and is not caused by changes of the cover layer, 
i.e., adsorbates4S We observed changes of the cover layer 
if we start with a base pressure of p = 10 mbar, but 
not for p < 10 _9 mbar. In Fig. Q] we plot cr°, a\, and a\ 
of YTi0 3 from 0.75 to 5.8 eV at 15 K. The data are con- 
sistent with the unpolarized room-temperature data of 
Ref. HO and with infrared transmittance and reflectivity 
results obtained in our groupi^l The latter revealed an 
onset of interband excitations at about 0.6 eV (see Fig. 
[2]). Recently, the effect of oxygen defects at the surface 
of YTi0 3 has been discussed^ We address this issue in 
the appendix. 



III. RESULTS 

Undoped YTi0 3 and SmTi0 3 are Mott-Hubbard in- 
sulators. In the ground state there is a single electron 
in the 3d shell at each Ti site. It is well accepted that 
the absorption above the gap corresponds to excitations 
from the LHB to the UHB, i.e., to the creation of an 
empty and a doubly occupied site, \d 4 d}) — > |d°<i 2 }. The 
strong increase of ai(u>) above «4.5eV (sec Fig. [IJ re- 



flects the onset of charge-transfer excitations from the 
2p band to the UHB, \d 4 p 6 } -> \d 2 p 5 ). The difference 
in spectral weight can be attributed to the Ti-0 hop- 
ping t p d- ci (w) oc t 2 d for charge-transfer excitations and 
(Ti(u>) oc t 4 d /A 2 for Mott-Hubbard excitations, where A 
denotes the charge-transfer energy. 

For YTi0 3 , photoemission and inverse photoemission 
spectroscop y 47 ' 53 ' 54 ' 55 yield A = 6cV and an on-site 
Coulomb interaction U ~ 5 eV^ where U denotes the 
Coulomb repulsion if both electrons occupy the same real 
orbital. In a single-band Hubbard model, the splitting 
between LHB and UHB is given by U (cf. inset of Fig. 
[T]). However, for a quantitative description of cr(ui) and 
for a reliable peak assignment one has to take all five 3d 
orbitals into account ^a^IL^ 



A. Multiplet assignment and Hubbard exciton 

Figures [3] and 2] focus on the inter-Hubbard-band exci- 
tations of YTi0 3 and SmTi0 3 below 4.5 eV. In YTi0 3 , 
three peaks are observed at 1.95 (A), 2.55 (B), and 
4.15 eV (C). In SmTi0 3 , we find two pronounced peaks at 
1.9 and 3.7 eV. Additionally, there is a shallow shoulder 
at 2.5 eV, particularly noticeable for the b axis. 

For a Mott-Hubbard insulator, one expects that a lo- 
cal multiplet calculation yields a reasonable assignment 
of the LHB-UHB excitations, 3 ^! 3 - 7 .^ In terms of local 
multiplets, the excited \d°d 2 } states can be distinguished 
according to the d 2 sector, because d° is an empty shell. 
The d 2 sector is split into a series of multiplets by the 
electron-electron interaction, the crystal field, and the 
hybridization with the ligands^ 7 - We start from cubic 
symmetry, in which case the crystal field and the hy- 
bridization give rise to a splitting of the 3c? orbitals into 
a triply degenerate t 2g level and a doubly degenerate 
e g level at higher energy. The splitting is denoted by 
10 Dq, which roughly can be estimated as 2 ± 0.5 eV 
[HHHliil. The electron-elect ron interaction within 
the 3d shell can be parameterized by the three Slater 
integrals F°, F 2 and F 4 . Values of ,F 2 =6.75eV and 
F 4 /F 2 w 5/8 are characteristic for d 2 Ti 2+ ions in a 
crystali 5 - 7 - The only parameter that can be adapted is F a , 
which drastically deviates in a solid from the ionic value 
due to screening effects. 

For F° = 3.60cV (or U w 4.5cV [H) the ^d 1 ) -> 
\d°d 2 } excitation energies are given in Fig. [51 focusing on 
the four multiplets lowest in energy: the triplet 3 Ti, the 
singlets 1 T 2 and 1 E, and the triplet 3 T2. For an intuitive 
picture we consider the strong crystal-field limit (10 Dq 
3> U), as sketched on the right hand side of Fig.[5j In this 
limit, there is one electron in the t 2g level and one in the 
e g level in the 3 T2 state, whereas both electrons occupy 
the t 2g level in the three other states. It is common 
to consider the simplified Kanamori scheme^ 7 - with the 
Hund on-site exchange coupling Jh = "^F 2 + ^^-F 4 , 
resulting in J H = 0.6 ± 0.1 eV for d 2 Ti 2 +. For U = 4 - 
5 eV, the Kanamori scheme predicts the lowest excitation 
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FIG. 3: (Color online) Optical conductivity of YTiOs below 
the onset of charge-transfer excitations, i.e., in the range of 
the lowest excitations from the lower to the upper Hubbard 
band. Peak B is attributed to the lowest multiplet, whereas 
peak A is identified as an excitonic resonance. Peak C reflects 
the lowest excitation to an e g orbital. 



into the 3 T\ triplet at U — 3Jh ~ 2 — 3eV, separated 
from the singlets x Ti and 1 E by 2Jh ~ 1.2 eV (reflecting 
Hund's rule) and from the 3 T 2 state by lODqra 2cV, 
in qualitative agreement with the result of the rigorous 
calculation shown in Fig. [5] 



1. YTiOs 

Figure [5] clearly shows that the 3 Ti state is the low- 
est multiplet, more than 1.2 eV below the next multiplet 
for any reasonable choice of 10 Dq. Thus the small split- 
ting of 0.6 eV between peaks A and B in YTiC>3 cannot 
be identified with the difference between the 3 T\ state 
and any other multiplet. We conclude that both peaks A 
and B arc related to excitations into the 3 7i state. Peak 
C can be attributed to the 3 T 2 state, since only excita- 
tions into triplet states are allowed from a fully polarized 
ferromagnetic ground state within an electric dipole ap- 
proximation. Excitations to the singlet states 1 T 2 and 
1 E require a spin flip and thus are suppressed, at least 
at low temperatures. 

In the following, we discuss three scenarios for the 
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FIG. 4: (Color online) Optical conductivity of SmTiOa 



splitting between peaks A and B: deviations from cu- 
bic symmetry, band structure effects, and an excitonic 
resonance. The deviation from cubic symmetry lifts the 
degeneracy of the t<x g orbitals and thereby also of the 3 Ti 
state. The t^g splitting was found to be «0.25eV in in- 
frared transmittance,— Raman scattering,— and RIXS 
measurements^ This is clearly too small to explain the 
splitting between peaks A and B.— 

Now we address the possible role of band struc- 
ture effects. Based on the actual crystal structure, a 
LDA+DMFT study of YTi0 3 by Pavarini et alM does 
not show a splitting of the lowest peak in a\{uj). For 
U—5 6V and Jff=0.64eV, this peak has been predicted 
at 3.3 eV, and the optical gap is expected roughly at 
1.5 eV. This large value of the gap suggests that a smaller 
value of U is more appropriate. Good agreement between 
the prediction for the lowest peak and the observed en- 
ergy of 2.55 eV of peak B can be obtained by assum- 
ing U « 4.3 eV. This value of U also yields a good de- 
scription of the optical gap. Moreover, it corroborates 
the validity of our local multiplet calculation discussed 
above (see Fig. EJ, which for [7=4.5 eV predicts the low- 
est peak at about 2.5 eV. We stress that it is unreason- 
able to identify peak A at 1.95 eV with the peak found in 
LDA+DMFT, since this would require to assume a still 
smaller value of U, resulting in a very small gap. Indeed 
the LDA+DMFT calculation finds a metallic state for 
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FIG. 5: (Color online) Left: Calculated energies for a Id 1 ^ 1 ) 
— » |d°ci 2 ) excitation with different d 2 final states in a cu- 
bic crystal-field.— The Slater integrals were chosen as F° 
= 3.60 eV, F 2 =6.75eV, and F 4 =4.55eV, corresponding to 
C/ = 4.5eV and J H « 0.6 eV. For fODq = the ionic mul- 
tiplet structure is obtained. For 10 Dq ~ 1.5 - 1.9 eV (grey 
area) the energy of 4.I5eV for peak C is well described by 
excitations into the 3 T2 state. Right: sketch of the orbital 
occupation in the strong crystal-field limit for the 3 Ti triplet 
(red), the l Ti and 1 E singlets (green), and the %i Ti triplet 
(blue). 



Also the LDA+DMFT study by Craco et alM for 
the ferromagnetic phase of YT1O3 finds a single peak 
in <Ji(u>). Based on the parameter values of £/ = 4.75eV 
and Jjj = 1.0 eV, Craco et al. attribute this peak to 
peak A observed at 1.95 eV in our data. However, Jh 
is not expected to deviate strongly from the ionic value 
of Jh = 0.6 ± 0.1 eV discussed above. The lowest peak 
in <ti(w) is located at about U — 3Jy, thus the choice of 
Jjj = 1.0eV strongly underestimates the peak frequency. 
We emphasize that both LDA+DMFT studies&S 4 - find 
a single peak in <7i(w). Both studies investigate an effec- 
tive Hamiltonian for the t2 g sector, i.e., excitations to the 
higher-lying 3 T2 multiplet (peak C) are not considered. 

Experimental data neither support a splitting due to 
band-structure effects. In photoemission (PES) data of 
YTiC>3 the LHB is a single peak w 1.3 eV below the Fermi 
level4!£^£5 In inverse PES on Yi^Ca^TiOa (x=0 [H 
and 0.4 - 0.8 (Hi) the UHB can be identified with the low- 
est peak or shoulder « 1.5 - 2eV above the Fermi level. 
Both PES and inverse PES agree with the LDA+DMFT 
result^ 3 - for U — 4 - 5eV. Finally, U w 5.3 eV has been 
derived from 2p core-level PES [H3| (see (5(| for the com- 
parison of parameters derived by different methods). 

Altogether, both theoretical and experimental results 
support our interpretation that the splitting of 0.6 eV 
between peaks A and B does not result from the band 
structure and that peak B at 2.55 eV is the dominant 
excitation. 

In contrast to (inverse) PES, the optical conductiv- 
ity reflects particle-hole excitations and thus is sensi- 



tive to interactions between the particle in the UHB 
(i.e. a double occupancy) and the hole in the LHB. 
These interactions are also neglected in the LDA+DMFT 
calculation s 6 ^ 64 of a\(u>). We therefore identify peak B 
at 2.55 eV as a particle-hole excitation in which the par- 
ticle and the hole are well separated, whereas peak A at 
1.95 eV is interpreted as an excitonic resonance, where 
the particle and the hole remain close to each other. 
Note that peak A does not lie below the gap, i.e., it 
is not a truly bound exciton, but a resonance within 
the continuum. As discussed above, a Hubbard exciton 
may arise due to the attractive Coulomb interaction be- 
tween the particle and the hole. The nearest-neighbor 
electron-electron repulsion V of the extended Hubbard 
mod.elS*ZiSi2i i 10 i 11 i 12 i 13 i 14 is equivalent to a particle-hole 
attraction — V [66[. We arc not aware of an accurate 
experimental value of V for the titanates, but it is rea- 
sonable to assume V < 1 eV. For the ID charge-transfer 
insulator SrCuC"2, a value of V « 0.6 eV has been de- 
rived from the comparison of the line shape of the exci- 
tonic resonance in <j\ (lo) with predictions from dynamical 
density-matrix rcnormalization group calculations for an 
effective extended Hubbard models More detailed theo- 
retical studies of the extended Hubbard model in 3D arc 
required to decide whether the nearest-neighbor Coulomb 
interaction is sufficient to explain the splitting of 0.6 eV 
observed between peaks A and B. 



2. SmTiOs 

The magnetic ground state of -RTiC"3 changes from 
ferromagnetic to antiferromagnetic as a function of the 
size of the R ions 1 71 1 72 This change is accompanied 
by a crossover of both the character of the distortions 
of the oxygen octahedra and of the orbital-ordering 
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The optical conductivity of the anti- 



ferromagnet SmTiOa is given in Fig. 01 focusing on the 
range of the Mott-Hubbard bands below the onset of 
charge-transfer excitations at about 4.5 eV. At 300 K we 
observe two pronounced peaks at 1.9 and 3.7 eV and a 
shallow shoulder at 2.5 eV, most evident for the b axis. 
The multiplet structure discussed above for YTiC-3 also 
applies to SmTi03. In particular, one does not expect 
appreciable changes of the Slater integrals, i.e., of the 
electronic parameters U and Jjj. This is corroborated by 
the LDA+DMFT calculations by Pavarini et al.^ which 
predict the same peak frequency for (J\{uj) in YTi03 and 
in the antiferromagnet LaTiC>3. Therefore we identify 
the shoulder at 2.5 eV in SmTi03 with peak B at 2.55 eV 
in YTi03, whereas the asymmetric peak at 1.8 eV is at- 
tributed to an excitonic resonance (see below), equivalent 
to peak A in YTi0 3 . 

In contrast to the peak energies, the band width or 
equivalently the hopping matrix elements are expected to 
change significantly, resulting from the different Ti-O-Ti 
bond angles. A larger band width of SmTiOs agrees with 
the observation from transmittance measurements^ 1 - that 
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FIG. 6: (Color online) Top: Leading edge of the optical con- 
ductivity ct\{uj) of SmTi03. Bottom: Temperature depen- 
dence of 0i(w = 1.85 eV) (crosses, left axis) and of the lead- 
ing edge uj e , which we define by a\{uj £ ) = 350(ficm) _1 (full 
symbols, right axis). Open symbols show Cb e corrected for the 
change of the spectral weight, i.e., a\{u) e ) = c ■ 350 (ficm) -1 
with c= o-S , (1.85eV,T)/cr 1 (1.85eV,60K). 



the gap in SmTiOa is about 0.2 eV lower than in YTiC>3. 
Additionally, an increase of the hopping amplitudes gives 
rise to an increase of the spectral weight. This is further 
enhanced by the change of the orbital ground state . 71 i 72 
Experimentally, we find an increase of TVeff from YTiO.3 
to SmTi0 3 of roughly 25%, 50%, and 100% for the a, b, 
and c axes, respectively. 

As discussed above for YTiC>3, an interpretation of 
the peak at 1.9 eV in terms of the lowest multiplet is 
hard to reconcile with the LDA+DMFT result^ unless 
excitonic effects are considered. An exciton interpreta- 
tion is supported by the temperature dependence of the 
peak frequency observed for the b and c axes, showing an 
anomalous softening with decreasing temperature and an 
anomaly at Tjy. For the b axis we demonstrate this soft- 
ening in Fig. [5] We focus on the frequency cu e of the 
leading edge, which we define as a\(uj e ) = 350(f2cm)~ 1 . 
This has the advantage that uj e and also its tempera- 
ture dependence can be determined more accurately than 
the peak frequency itself. The disadvantage is that a 
softening of lu c in principle can be caused not only by 
a softening of the peak frequency, but also by an in- 
crease of either the spectral weight or the line width, 
and by a change of the line shape. We find a jump-like 
decrease of uo e at Tn, see lower panel of Fig. O This 
cannot be attributed to an increased line width, since 
the thermal contribution to the line width is expected 
to decrease with decreasing temperature. Moreover, for 
spin- carrying particles one expects that the band width 
is reduced upon entering the AF ordered state, thus the 



gap is expected to harden. For an estimate of the T de- 
pendence of the spectral weight we consider the value 
of <7i(lu) at the peak frequency. We find an increase 
of crJ(tt>=1.85eV) upon cooling below Tjv (see crosses in 
bottom panel of Fig. [6]), indicating an increase of the 
spectral weight. We use this T dependence of the ab- 
solute value to determine a corrected frequency of the 



leading edge, u/ e , defined as <Tj(w e ) 



350(f2cm)- 



with c = crJ(1.85eV,T)/crJ(1.85eV,60K) (open symbols 
in Fig. 0. This shows that the shift of w e is not due 
to the change of the spectral weight, it is caused mainly 
by a softening of the peak frequency or a change of the 
line shape. Both can be rationalized by the attractive in- 
teractions responsible for exciton formation, pulling the 
spectral weight to lower frequencies. 

The excitonic binding due to the kinetic energy is en- 
hanced in the AF ordered state, as discussed for 2D com- 
pounds in the introduction. Remarkably, the peak fre- 
quency of 1.95 eV is independent of temperature in the 
ferromagnet YTiC>3. It is an interesting question whether 
this difference between the two compounds arises from a 
change of the screening of nearest-neighbor Coulomb in- 
teractions, of the magnetic ground state, or of the orbital 
ordering pattern. A decisive identification of the driv- 
ing force for exciton formation requires further theoreti- 
cal investigations of the extended multi-orbital Hubbard 
model in 3D. 

The peak at 3.7 eV coincides with the minimum of 
(7i (w) observed in YTi0. 3 . This peak can be attributed 
to the lowest singlet multiplet ( 1 T2 and 1 E in cubic sym- 
metry). Due to the spin selection rule, the excitation to 
the singlet state is suppressed in ferromagnetic YTiC>3, 
but it is allowed in antiferromagnetic SmTiOs . This fea- 
ture is expected at about 2Jh ~ 1.2 - 1.3 eV above the 
lowest triplet peak, providing further support for the as- 
signment of peak B at 2.5 eV and the excitonic character 
of the peaks at 1.8 eV in SmTi0 3 and 1.95 eV in YTi0 3 . 



B. Hubbard exciton and orbital order 

For a 2D Mott-Hubbard insulator with AF exchange J 
on a square lattice, exciton formation is governed by the 
kinetic energyi 18 ' 19 i 20 i 21 i 22 The motion of a single parti- 
cle is hindered by the interaction of its spin with the AF 
background. This can be described in terms of a spin po- 
laron. Hopping of the bare particle on the energy scale 
t results in a trace of misaligned spins. Coherent motion 
of the dressed polaronic quasiparticle requires the emis- 
sion of magnons, i.e., the bare band width ~ t is reduced 
to the polaronic band width ~ J, which corresponds to 
an increase of kinetic energy. In this case, the kinetic 
energy is lowered by the formation of spinless excitons, 
which recover a larger band width. 

This mechanism may contribute to exciton binding 
in antiferromagnetic SmTiOs, but not in ferromagnetic 
YTi03. It is promising to investigate whether a simi- 
lar mechanism is at work in the case of antiferro- orbital 
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the hole heals out the trace of excited orbitals, see Fig. 
[7^ and f. Therefore, the motion of the exciton is not hin- 
dered by the antiferro-orbital order, and the exciton can 
hop on a larger energy scale than the hole or the double 
occupancy individually. Thus exciton formation here is 
equivalent to a gain of kinetic energy. 

More detailed knowledge on the value of the nearest- 
neighbor Coulomb interaction V and its relationship to 
the binding energy in 3D Mott-Hubbard insulators is 
required to decide whether this mechanism is realized 
in YTiC>3. The orbital ordering pattern in YTiC>3 is 
more complex than simple antiferro-orbital order ; 32 ' 33 
and hopping between orbitals of different type is not ex- 
actly zero. Still Fig. [7J may be relevant for the ab plane, 
since hopping from the lowest orbital on one site (a "cir- 
cle" in Fig. [7J to the lowest orbital on a neighboring site 
(equivalent to a "square" ) is 2-3 times smaller than hop- 
ping to the excited states! 60 ' 63 



FIG. 7: Sketch of the suggested formation and propagation 
of a Hubbard exciton (dashed line). We consider two types of 
orbitals (circles and squares, e.g., d xy and d xz ) per site, where 
hopping is zero between orbitals of different type (crossed 
out arrows). Full (open) symbols denote occupied (empty) 
orbitals. (a) Ground state with antiferro-orbital order, (b) 
Creation of a hole and a double occupancy on sites 2 and 
3, respectively, (c)-(f) Propagation of the double occupancy, 
the hole, or of an exciton (see main text for more details). 



order. For illustration and simplicity, we consider a ID 
model with two orbitals per site, e.g., d xy and d xz for a 
chain running along the x direction. In Fig. [JJ the two 
types of orbitals are denoted by circles and squares, re- 
spectively. Hopping between neighboring sites is allowed 
only between orbitals of the same type; it is zero be- 
tween orbitals of different type. Black and grey symbols 
in Fig. [JJ refer to occupied orbitals, whereas empty sym- 
bols denote empty orbitals. The ground state in Fig. [7^, 
exhibits antiferro-orbital order, i.e., xy (circles) and xz 
orbitals (squares) are occupied in an alternate fashion. 
The empty orbitals are at higher energies due to, e.g., 
the ligand-ficld splitting. An excitation from the LHB to 
the UHB, i.e., ^d 1 ) -> \d n d 2 ), is illustrated in Fig. [7p. 
Site 2 is empty, and site 3 is doubly occupied. The mo- 
tion of the double occupancy to sites 4 and 5 is depicted 
in Figs. [7J; and d, respectively. The central point is that 
this motion leaves a trace of orbitally excited states, i.e., 
on sites 3 and 4 the energetically unfavorable orbitals are 
occupied (grey symbols). This results from the restric- 
tion that hopping is only allowed within the same type 
of orbital. As discussed above for the case of spins, co- 
herent motion of the quasi-particle requires the emission 
of orbital excitations, in our example the de-excitation of 
sites 3 and 4. Therefore, the band width is reduced from 
the bare band width ~ t to the energy scale of the or- 
bital excitations, corresponding to an increase of kinetic 
energy. However, if the hole accompanies the double oc- 
cupancy forming an exciton (dashed line), the motion of 



C. Temperature dependence of the spectral weight: 
spin and orbital selection rules 

The spectral weight is determined by the spin and or- 
bital selection IU \ es 4^6^8^9j^^3^^ There . 

fore, the T dependence of the spectral weight is expected 
to reflect changes of the spin-spin correlations and/or 
of the orbital occupation. By considering the nearest- 
neighbor spin-spin correlations, the absolute value of 
the spectral weight has been calculated for instance for 
LaMnOa HHI^l and LaSrMn0 4 H^. These calcula- 
tions yield a convincing description of the experimental 
results, the maximum difference is less than a factor of 2. 
For a 3D magnet one expects that the spin-spin correla- 
tions are small above the ordering temperature. In fact, 
the change of the spectral weight above T/v is small in 
the 3D antiferromagnet LaMn03^ In contrast, the 2D 
antifcrromagnet LaSrMn04 with a Neel temperature of 
Tyv = 130 K, exhibits a significant T dependence of the 
spectral weight up to 300 K, which can be attributed to 
enhanced quantum fluctuations in ZD 2^ 

For the lowest excited triplet state ( 3 Ti in cubic no- 
tation) of orbitally ordered YTi03, Oles et al£^ pre- 
dicted a change of 25% of the spectral weight between 
the paramagnetic and the ferromagnetic state in the ab 
and c directions. This can be understood by the evolu- 
tion of the nearest-neighbor spin-spin correlation func- 
tion (Si ■ Sj + 3/4), which equals 1 in the ferromagnetic 
state and 3/4 in the paramagnetic state, i.e., a redistri- 
bution of 25%. 

We analyze the integrated spectral weight in terms of 
the effective carrier concentration N e g, 

n cs =^ r\^)^ (i) 

where tu c i and tu C 2 denote the frequency range of interest, 
m is the free electron mass, e the elementary charge, and 
1/Vq the density of Ti ions. For Vq we use the value 
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FIG. 8: Temperature dependence of the effective carrier con- 
centration iVefi (see Eq. [T} of YTiOa for uj c i — 1.6 eV and 
uj C 2 = 2.6 eV (left) and for lj c i = 2.6 eV and a> C 2 = 3.9 eV 
(right). 
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observed at 290 K, which differs from the 2 K value by less 
than 1%M The T dependence of N eS is given in Fig. [S] for 
YTi0 3 and in Fig.[9]for SmTi0 3 . In YTi0 3 , the spectral 
weight increases between 1.6 and 2.6 eV (left panel of 
Fig. [H]) in the a and b directions upon cooling down from 
room temperature. We find an anomaly in the vicinity 
of T c , i.e., an additional increase of spectral weight with 
decreasing temperature. This additional increase starts 
at about 1.5 - 2 T c and amounts to less than 5% below 
50 K, much smaller than predicted. At the same time, 
one finds an anomalous decrease of spectral weight with 
decreasing temperature between 2.6 and 3.9 eV, again in 
a and b. Moreover, the spin selection rule cannot explain 
that the spectral weight in this 3D ferromagnet shows a 
strong T dependence up to 300 K > 10 T c . 

In SmTiOa, we also find pronounced anomalies in the 
vicinity of the magnetic ordering temperature, Tjy = 
53 K. The strongest change of N e g of up to 10% is ob- 
served in the c axis, significantly larger than in YTi0 3 . 
In SmTi0 3 , the sign of the change at T/v depends on 
the polarization (sec left panel of Fig.O, which certainly 



cannot be attributed to the spin selection rule in a G- 
type antiferromagnet. In comparison to the excellent 
agreement found between experiment and theory in the 
manganiteS ) 35 ' 38 this failure appears as a puzzle. 

Alternatively, we consider the orbital selection rule. 
In the manganites, the orbital occupation has been as- 
sumed to be independent of T due to the large e g split- 
ting of roughly 1 eV^^ In both YTiOa and SmTi0 3 , 
the t 2g splitting is only w 0.25 eV (mUlEl, 

opening the 

possibility for small changes of the orbital occupation 
as a function of T. A change of the orbital occupa- 
tion affects the effective Ti-Ti hopping amplitude i, with 
A^off (xt 2 (x tp d /A 2 . An increase of the occupation of the 
planar xy orbital may for instance give rise to an increase 
of the spectral weight within the xy plane accompanied 
by a decrease of spectral weight along z. Thus a change 
of the orbital occupation at TV can very well account 
for the observed polarization dependence. Based on a 
detailed analysis of the crystal structure, thermal expan- 
sion and magnetostriction of i?Ti0 3 , Komarek et alM- 
conclude that magnetism affects the crystal structure, 
which in turn drives a change of the orbital occupation. 
Remarkably, the shape of the oxygen octahedra changes 
significantly as a function of temperature, whereas the 
variation of the tilt and rotation angles is small* 3 ^ Both 
the lattice distortions and the orbital occupation adapt 
in order to enhance the gain of energy within the spin 
system. The effect is most pronounced at the magnetic 
ordering temperature, but extends also to higher temper- 
atures, in agreement with our data. Moreover, Komarek 
et alM: pointed out that the change of the orbital occupa- 
tion is significantly stronger in SmTi0 3 than in YTi0 3 , 
again in agreement with our results. The occurrence of 
pronounced effects in SmTi0 3 is attributed to the fact 
that SmTi0 3 is close to the crossover from antiferromag- 
netic to ferromagnetic order^i In the optical data, the 
effect of the orbital selection rule possibly overrules that 
of the spin selection rule, which appears as a failure of 
the latter. 

Note that the change of the lattice constants at T/v can 
only account for a change of N e s of the order of 1%. This 
estimate is based on the Harrison rules^ for the hopping 
amplitudes. We emphasize that binding phenomena such 
as the formation of excitons or resonances in general are 
very sensitive to temperature. With decreasing temper- 
ature, the attractive interactions responsible for the ex- 
citon formation pull down the spectral weight to lower 
energies, in agreement with the change of the line shape 
observed in YTi0 3 (see Fig. [3]) and the shift of the ab- 
sorption edge of SmTi0 3 discussed above (see Fig. [5]). 



D. Anisotropy 

In order to understand the anisotropy observed in 
YTi0 3 between the ab plane and the c direction, we 
address the matrix elements for the optical excitation 
| cZ 1 cZ 1 ) — > \d°d 2 }. Our Hamiltonian includes the crys- 
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tal field, the on-site Coulomb correlations of the d 2 con- 
figurations, and the hopping between the two Ti sites 
(for details, see Ref. 160 ). We find realistic values of 
the exchange coupling constants for the different direc- 
tions as well as an orbital ground state which is in excel- 
lent agreement with x-ray, neutron and other theoretical 
resultsi 32 ' 33 ! 51 ' 63 ! 67 From the effective Ti-Ti hopping ma- 
trices t ab and t c one can estimate the anisotropy of the 
spectral weight from 

= N? s + N b cS = v (f#) 2 + (ffi) 2 

where denotes the effective hopping matrix clement 
between the t^g orbitals i and j on adjacent sites. We 
obtain N^g/N^ s « 5.1. Using the hopping matrices pub- 
lished by other groups, we find a value of 1.1 [68| or 
3.5 (63|. In the next step the optical conductivity has 
been calculated using the Kubo formula. We assume a 
fully polarized ferromagnetic ground state. We address 
only excitations into the lowest triplet state with a t\ g 
configuration, because here the point-charge approxima- 
tion gives reliable results. We predict N£g/N£ s w 4.5. 
The small difference to the value of 5.1 derived from 
the simplified approach considered in Eq. [2] arises be- 
cause here the energies of the excited states and the Ti- 
Ti distance are taken into account. For lu c i = 1.6 eV 
and uj C 2 = 2.6 eV (see Eq. [T]), we experimentally find 
N^/N* s « 2 (see Fig.©, within the range predicted by 
the different theoretical approaches. 



is much smaller than predicted, whereas in SmTi03 even 
the sign of the temperature dependence disagrees for cer- 
tain polarization directions, a puzzling result. However, 
a small change of the orbital occupation at the magnetic 
ordering temperature^ can account for the polarization 
dependence and also explains the larger temperature de- 
pendence found for SmTi03. In contrast to the man- 
ganites, such a change of the orbital occupation is feasi- 
ble in i?TiC>3 because the t2 g splitting amounts to only 
0.25 eV. Furthermore, the increase of spectral weight at 
low frequencies with decreasing temperature is in agree- 
ment with an exciton scenario, since binding phenom- 
ena are expected to exhibit a strong temperature depen- 
dence. An increased binding at low temperatures pulls 
down spectral weight to lower frequencies and also ex- 
plains the anomalous softening of the loading absorption 
edge observed in SmTi03. 

The importance of excitonic effects for the description 
of (t(lo) is well established for low-dimensional correlated 
insulators. The attractive interaction responsible for ex- 
citon formation arises from a gain of either Coulomb or 
kinetic energy. We have pointed out that exciton forma- 
tion may lower the kinetic energy in an orbitally ordered 
state. Our results call for further theoretical studies of 
exciton formation in the extended multi-orbital Hubbard 
model in 3D. A quantitative description of this binding 
phenomenon is essential for a consistent explanation of 
optical and photoemission data and will provide impor- 
tant information on electronic correlations in Hubbard 
systems. 



IV. SUMMARY AND CONCLUSIONS 

In summary, we report on optical excitations from the 
lower to the upper Hubbard band in the ferromagnct 
YTiC>3 and in the antiferromagnet SmTi03. At 15 K 
we find peaks in the optical conductivity <j\(ui) at 1.95, 
2.55, and 4.15cV in YTi0 3 and at 1.8 and 3.7cV in 
SmTiOs, which also exhibits a shallow shoulder at 2.5 eV. 
For these Mott-Hubbard insulators, a local multiplet sce- 
nario is expected to yield a reasonable peak assignment, 
as reported for the manganitesi 35 i 38 For ?7«4.5eV and 
Jh = 0.6 ± 0.1 eV, our local multiplet calculation offers 
a quantitative description of the peak positions at 2.5, 
3.7 and 4.15 eV. The peak at about 2.5 eV is attributed 
to excitations into the lowest d 2 multiplet ( 3 Ti in cubic 
symmetry) with an energy of roughly U — 3 Jh [62J- The 
peak at 3.7 eV corresponds to the lowest d 2 singlet states, 
and the peak at 4.15 eV is attributed to the lowest state 
with a t\gGg configuration. This assignment is in agree- 
ment with photoemission and LDA+DMFT results. The 
peaks at 1.95 eV in YTi0 3 and 1.8 eV in SmTi0 3 are 
interpreted in terms of an excitonic resonance, thereby 
explaining their low energy. 

The temperature dependence of the spectral weight 
disagrees with predictions based on the spin selection 
rule. In YTi03 the observed temperature dependence 
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Appendix: Role of oxygen defects 

Recently, Kovaleva et al£^ studied YTi03 by ellipsom- 
etry and reported on complications which they attribute 
to oxygen defects arising from polar surfaces. They ob- 
served peaks in a\(u>) at 1.95, 2.9, and 3.7 eV. The over- 
all temperature dependence observed in Ref. [13 is very 
weak, showing a crossover at 100 K, but no anomaly at 
T c within the experimental accuracy. In the frequency 
range studied by us, the main effects of oxygen defects 
were identified as (i) a shift of the fundamental absorp- 
tion edge to lower frequencies, (ii) an absorption band at 
about 0.8 eV, (iii) the absence of a pronounced minimum 
at 4.5 eV, and (iv) a shift of the onset of charge-transfer 
excitations to lower frequencies. These shifts have been 
attributed to localized states at the edge of the electronic 
bands. This sensitivity of the fundamental absorption 



10 




Energy (eV) 

FIG. 10: (Color online) Comparison of the optical conductiv- 
ity for the b axis of YTi03 with spectra reported by Kovaleva 
et al.r— which were determined from either the ab or the be 
surface of a freshly polished sample, or from the ab surface of 
a sample measured after 1.5 years. 



edge to doping away from the half-filled Mott insulator 
has been studied in Y^Ca^TiOg \jx$7(%. In Fig. [TO] 
we compare our data with the results of Ref. Our 
data show both the largest fundamental absorption edge 
and the largest onset frequency for charge-transfer ex- 
citations, in combination with a pronounced minimum 
at 4.5 eV. We find good agreement with the spectrum 
of &i(uj) below the fundamental gap determined in our 
group by transmittance and reflectance measurements 
on thin single crystals^ (see Fig. [2]). The transmit- 
tance clearly reveals bulk properties. These data show no 
defect-induced absorption below the gap, the single weak 
feature observed at about 0.3 eV has been undoubtedly 
identified as a phonon-activated orbital excitation . 51 ' 59 ! 61 
Finally, we find clear anomalies in the vicinity of the 
magnetic ordering temperatures, both in YTiC>3 and in 
SmTiC>3 (see Figs. [8] and [9]). The combination of all these 
observations provides strong evidence that we have ob- 
served the intrinsic properties of YTiC>3. 
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